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In this publication we investigate dynamics of a flat FRW cosmological model with a non- 
minimally coupled scalar field with the coupling term ^Rip'^ in the scalar field action. The 
quadratic potential function V{ip) oc ip'^ is assumed. All the evolutional paths are visualized 
and classified in the phase plane, at which the parameter of non-minimal coupling ^ plays the 
role of a control parameter. The fragility of global dynamics with respect to changes of the 
coupling constant is studied in details. We find that the future big rip singularity appearing 
in the phantom scalar field cosmological models can be avoided due to non-minimal coupling 
constant effects. We have shown the existence of a finite scale factor singular point (future 
or past) where the Hubble function as well as its first cosmological time derivative diverge. 
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INTRODUCTION 



Recently scalar fields have played a very important role in cosmology. They are used in many 
phenomenological models like quintessence jl|,y]- Scalar fields are also very important in description 
of dynamics in the loop quantum cosmology, which base on the background independent theory 
without the canonical notion of time. In this theory one scalar field is chosen as an internal clock 



for other fields 



The scalar fields with a potential function are also very important in modelling 



of inflation. For example a scalar field with the simplest quadratic potential function was assumed 
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in Linde's conception of chaotic inflation [4]. The 5 years of WMAP observations rejected many 
inflationary scenarios (potential functions V{tl^)), while models with a simple quadratic potential 
are admitted on the la confidence level. 

In the standard quintessence energy density of the minimally coupled to gravity scalar field 
mimics the effective cosmological constant. Of course the detailed evolution is dependent upon a 
specific form of the potential V{'ip) but the V'^ contribution can be treated as a leading order term 
of expansion of the potential function. 

We extend the quintessence scenario by incorporating the non-minimal coupling constant 
Q, y, y, [3]. In this paper we present a phase space analysis of the evolution of a spatially flat 
Friedmann-Robertson- Walker (FRW) universe containing a non-minimally coupled to gravity scalar 
field, both canonical and phantom, with the simplest form of quadratic potential function. The 
similar analysis for the case of minimally coupled scalar field was performed in We extend this 
analysis on models with a non-zero coupling constant ^ which plays the role of a control parameter 
for an autonomous dynamical system on the phase plane. Therefore the location of fixed points 
(physically representing asymptotic states of the system) as well as their character depends upon 
the value of The values of parameter ^ for which the global dynamics changes dramatically are 
called bifurcation values. In our previous paper [l^, in case of conformal coupling and quadratic 
potential function we have shown that phantom cosmology can be treated as a simple model with 
a scattering of trajectories whose character depends crucially on the sign of the potential function. 
We also demonstrated that there is a possibility of chaotic behavior in the flat Universe with a 
conformally coupled phantom fleld in the system considered on the non-zero energy level. 

The minimally coupled scalar field endowed with a quadratic potential function has a strong 
motivation in infiationary models and its generalizations with a simple non-minimal coupling term 
have been studied {l3] in the context of the origin of the canonical inflaton fleld itself. The 
physical motivation to investigate a non-minimally coupled scalar field cosmological models could 



be possible application of this models to inflationary cosmology or to t 



for example 
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le present dark energy (see 



The coupling constant ^ is a free parameter of the theory which should be estimated from the 
observational data. Recently we have shown that distant supernovae can be used to estimate the 
value of this parameter [l^ . 

The main advantage of dynamical system analysis is that we can visualize all the trajectories 
of the system admissible for all initial conditions. Therefore one can classify generic routes to the 
accelerating phase (the de Sitter attractor where = —p^)- This attractor corresponds to the 
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model with the cosmological constant. 

The paper has a following organization: in section [TI] we reduce dynamics to the form of an 
autonomous dynamical system which describes both canonical and phantom scalar field models. 
Section IIIII is devoted to a detailed analysis of the phase portraits for different values of the 
parameter ^. In this section we also discuss the change of evolutionary scenarios upon the value of 
parameter ^. 

II. NON-MINIMALLY COUPLED SCALAR FIELD COSMOLOGIES AS A 

DYNAMICAL SYSTEM 

We assume the flat model with the FRW geometry, i.e. the line element has the form 

ds^ = -df + a^it)[dr^ + r^{de^ + sin^ 9d(p'^)], (1) 

where < (f < 27r, < 9 < tt and < r < oo are the comoving coordinates and t stands for the 
cosmological time. 

It is also assumed that a source of gravity is a scalar field ip with a generic coupling to gravity. 
The gravitational dynamics is described by the standard Einstein-Hilbert action 

Sg = J d'xV^R, (2) 

the action for the matter source is 

^^ = -11 d'x^ [eig^'i^^i^u + ^R^^) + 21/(V')] ■ (3) 

where rrip^ = l/(87rG) = 1/ k and R = 6(a/a + a? /a?) and e = +1, —1 corresponds to the scalar 
field and the phantom scalar field, respectively. For simplicity and without lost of generality we 
will assume 47rG/3 = 1 which corresponds to k = 6. 

After dropping the full derivatives with respect to time we obtain the dynamical equation for 
scalar field from variation 8{Sg + S^)/5ip = 

^ + 3Hij + ^R^P + eV'{'tp) = 0. (4) 

as well as the energy conservation condition from variation 5{Sg + S^)/5g = 

If we include other forms of matter this condition can be expressed as 

3 

-H'^ = Pr + Pm (6) 
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where pr and pm are energy densities of radiation and matter, respectively. It can be shown that 
for any value of ^ scalar field behaves like a perfect fluid with energy density p^ and pressure 



Pi> = - ^(V') - [2^(V'')- + (V^'J^ (2i? + 2,H^)^^] . (7b) 



Changing the dynamical variables according to the relation 

■ dV' d dip , 
11) = — = -— — = HiIj 
at a d m a 

we can express the Hubble function as 

^2 ^ 2 ^(V') +Pr+Pm 

f -e[(l-60V''2 + 6e(V'' + V')2]' 
The denominator of ([8]) equal zero denotes a line in the phase space of singularities of the Hubble 
function which separates the phase space in two regions one physical > 0, and the second one 
nonphysical < 0. It does not depend on the form of the potential function but only on a value 
of the coupling constant. 

The Euler-Lagrange equations for the system under consideration are given in the form 



TT + = ~2a— — 

d?7^ dry^ dr] drj 



a^-T-T + eCaV'TT = -"^a—-- - ea*F'(V'), (9a) 



, 2(1 - ^'^Cr ) - ed^ai;-^ = -ea{l - 6^ -f + el2^^—-f- + 4a^y(V') + Pm,o- (9b) 
K dr/^ d??^ V dry / d?? dry 

where stands for the conformal time, dr] = dt/a. 

After the elimination of the scale factor and its derivative from system ([9]) we obtain the 

condition 

'6 



(< + -0') (- - £6^(1 - 60V'^) -eiP''^{l-Qi){iP' + Qi^p)- 



(10) 



1 r 6 

-< F — 

where the prime denotes differentiation with respect to the natural logarithm of the scale factor 
Now we can simply present equation pop in the form of the autonomous dynamical system 



V-' = y, (lla) 



^(1 _ 60 (y + eev-) 1 (f - 4(1 - 6^2/' + 6e(y + ^) 

y' = -y+ " 



X (e^y'(^)(i-eKCV'(y + V')) + (4F(V')+p^)(y + 6eV'))- (nb) 
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TABLE I: Finite critical points and their characters. 



^0 = (I) e6C(l-64) ' = ~^T^15^ 

'^O = (f ) £6;(l- ~60 ' 

yo = -4eY34|'/'0 ± ^ 3(1-44) 



existence 

-1 e = 



-1 



< ^ > 

0<e<H<Oore>i 

6 ^ "5 ^ 3 



eigenvalues 



Ai 



(!) (^3-5)^, A. = i(|)^(e3 + 5)^ 



^(1-30(3-250 
Ai = (f )2/o, A2 = (f )2yo 
Ai--(|)2yo, A2 = (f)6^V^o 
Ai = -(f)2yo,A2 = (f)2yo 



In what follows we will assume a quadratic potential function V{il)) = l/2m^'0^) ^-iid that there is 
no other form of matter than the scalar field, i.e that pr = Pm = 0. It is easy to notice that in 
this case the dynamics does not depend on the change of a sign of the potential F(V') — y(V')- 
Finally, the dynamical system is in the form 



(12a) 



a 



dy 

da 
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yi'(-- £6^(1 - GOr) + e(i - 60W(y + ecV')- 

- - 4(1 - 60y^ + 6e(y + V')'] ) - SK^^iy + ^P)) + 2V'(y + G^i^)) . (12b) 

where we have made the following "time" transformation 



where the parameter a 



a± = ^p^^- £6^1 - t/tT 
dfj \n /dm (a) 



a = +1 V(- - £6^(1 - 60-0^) > 0, 

a = -1 ^ V(- - £6^(1 - 60V'^) < 0. 



(13) 



(14a) 
(14b) 



was introduced to preserve the orientation of the trajectories in this way that on all of the phase 
portraits direction of arrows indicate the direction of growth of the scale factor. 



III. PHASE SPACE ANALYSIS OF DYNAMICS 



In this section we present detailed discussion of the type of critical points of the dynamical 
system p^ . 
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a) 



b) 




FIG. 1: The phase portraits for the canonical scalar field e — 1 for: a) minimal ^ = and b) conformal 
^ = 1/6 coupling. The shaded region denotes nonphysical part of the phase space for the strictly positive 
potential function. If the potential function is strictly negative the meaning of the regions is reversed. The 
shape of the border between the regions does not depend on the shape of the potential function. At the 
border of the physical region we have two symmetric critical points at the ip axis for both cases. The value 
of at that points is finite. The presence of a saddle type critical point in the case b) is the effect of 
non-zero ^. 

If we write the evolutional equations for the non-minimally coupled scalar field in the form of 
the dynamical system, the first step would be an identification of the critical points of the system. 
Physically they represent asymptotic (or stationary) states of the system under considerations and 
mathematically correspond to vanishing r.h.s. of the system. The second step is a characterization 
of the type of critical point which can be performed after calculation of the eigenvalues of the 
linearization matrix calculated of this critical point. The critical points are usually represented by 
physically interesting solutions and these solutions can be attractors for trajectories in its neigh- 
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FIG. 2: The phase portrait for the canonical scalar field e — 1 and coupling constant < ^ < 1/6. The 
shaded region is nonphysical: < for V{ijj) > 0. There are three types of critical points at the finite 
domain of the phase space: 1) i/iq = 0, j/q = 1 and ~ const, which are of a stable or unstable node type; 
2) ipQ = 1/6^, Uo = 0, = oo of a saddle type 3) ipo 0, yo ^ 0, = oo of unstable node type for 
F(V') > and stable node type for V{ip) < (shaded region). The dashed line denotes singularity of "time" 
transformation (|13p . In comparison with the phase portrait from Fig. [T] for conformal coupling we can see 
that both phase portraits are equivalent at the physical domain. 

borhood which evolve to it independently on the initial conditions. In the quintessence cosmology 
we are looking for the attractors, which give rise to solutions with desired property but we would 
like to know whether it is typical (generic) solution or exceptional (non-generic). This is a reason 
of our interest in the stability of the critical points. 

For full dynamical analysis investigation of the behavior of the trajectories at infinity is needed. 



30|. If we complete the phase plane by 



It can be performed by construction of Poincare sphere 
a circle at infinity which is a projection on the equator, then we receive the global phase portrait 
with a circle at infinity. In our case r.h.s. contain the non-minimal coupling constant ^ as a free 
parameter. The global phase portraits depend on the value of this parameter but for some ranges 
of the values of phase portraits can be equivalent (indistinguishable from the dynamical system 
point of view). If we fix the value of non-minimal coupling, then one can study the influence of 
this parameter on the global dynamics. The equivalence of the phase portraits is established by 
homeomorphism preserving direction of time along the trajectories. 

Critical point: tpo = 0, Vq = £ exists only for the canonical scalar field e = +1. Direct calculation 
of the Hubble function ([8]) at this point gives an undefined symbol It is why we use the linearized 
solutions in the vicinity of this critical point to show that the Hubble function at this point is finite 
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a) 



b) 



c) 




FIG. 3: The phase portraits for the canonical scalar field e = 1 and for the specific values of coupling 
constant: a) ^ = 3/16, b) ^ = 1/4, c) ^ = 3/10. In the cases a) and b) there are the critical points at infinity 
of a mixed type (multiple critical points). At the physical domain the phase portraits are topologically 
equivalent. 



9 



a) 



b) 



c) 




FIG. 4: The phase portraits for the canonical scalar field £ = I and values of coupling constant: a) 
1/6 < ^ < 3/16, b) 3/16 <S,< 1/4, c) C = 1/3. 
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FIG. 5: The phase portraits for the canonical scalar field e = 1 and negative coupling constant ^ < 0: in 
figure b) ^ is greater, but still negative, than in figure a). In the limit ^ — > 0^ we receive phase portrait in 
Fig. [TJi. 



and depends on the initial conditions of the hnearized solutions. They are in the form 

il>{a) = ^/^(i) exp (aAiO"), (15a) 
y(o-) = yo- eCV'Ci) exp {aXia) + (e^V'Ci) + (2/(i) - yo)) exp (aA2cr), (15b) 

where Ai = yo and A2 = 2yo are eigenvalues of the linearization matrix calculated at this critical 
point, ^/'(i) and are initial conditions and yo is a coordinate of the critical point. 

Inserting those solution into the formula dS]) we receive that the Hubble function in the vicinity 
of the critical point (ipQ = 0, i/q = e) is 



Hiin = rn^i'li) exp (2aAicr) - 6^(1 - e^V^i) exp {2a\ia) - 22/0(6^^^(1) + (^{i) - yo)) exp (QA2cr)- 



(36CV^i) + 12CV'{i)(y{j) - yo) + {y(i) - yof) exp (2QA2cr) 



(16) 
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FIG. 6: The phase portraits for the phantom scalar field e = — 1 and: a) minimal ^ = and b) conformal 
f = 1/6 coupling. All the phase space is admissible only for the positive potential function. We can 

conclude, that for negative potential functions in the case of minimally or conformally coupled phantom 
scalar fields, the scale factor is not a monotonic function of the cosmological time. For the case b) a global 
attractor represents the de Sitter state with w^, = — 1. There arc two types of trajectories which tend to this 
attractor: 1) trajectories starting from ip — 0, tp' ioo state, and 2) two single trajectories representing a 
separatrix of saddles at infinity (not shown). 



then taking the limit value of this function for a — > ±00 (depending on the critical point yo = =Fl) 
we receive 



limifii = m2- 



(i) 



(17) 



-2yo(6eV(.) + - yo)) - 6^(1 - QOi^fi) 
which is always a positive quantity. For the special values of minimal ^ = (see Fig. [1^) and 
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FIG. 7: The global phase portraits for the phantom scalar field and values of coupling constant: a) < ^ < 
3/25 and b) 3/25 < < 1/6. In the case a) in the finite domain the critical domain is of a stable node type 
and in the case b) of a focus type. 



conformal coupling ^ = 1/6 (see Fig. [Tb) the values of the Hubble function are 



-2yo(V'(i) + {y{i) - yo)) 



for ^ = 0, 



for f = -. 

^ 6 



Critical point: iJjq 



1 

£6^(1-60 



yo 



-e 



1 



(i-eov-o 



is very interesting because the Hubble function 



at this point is singular = oo and H = (hH^)' = oo. 



13 




FIG. 8: The global phase portraits for the phantom scalar field e ~ —1 and for the specific values of coupling 
constant: a) ^ = 3/16, b) f = 1/4, c) ^ = 3/10. In cases a) and b) one of the critical points at infinity is of 
a mixed type (multiple critical points) (see Fig. [3]). On all figures one can see trajectories starting from the 
unstable node and landing at the stable focus as a generic scenario of route to the de Sitter state. 
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FIG. 9: The global phase portrait for the phantom scalar field e = — 1 and distinguished value of coupling 
constant ^ — 1/3. In this case critical point at the finite domain of the phase space is located at the line of 
singularities of the time transformation . 




FIG. 10: The global phase portrait for the phantom scalar field e = — 1 and the values of coupling constant 
> 1/3. The characteristic critical point of a focus type disappeared. 

Linearized solutions in the vicinity of this critical points are 

ijj{a) = V'o + exp (aA2(T) - ipo) , (18a) 
yicr) = yo + exp (aAifj) (^2(1 - 3^)(V'(i) - V'o) + (vii) - yo)) - exp (aA2a) (^2(1 - 3^)(V'(i) - V'o)) 

(18b) 

where Ai = — yo ^md A2 = — 2yo are eigenvalues of the linearization matrix at the critical point, 
'(/'(j) and are initial conditions and ipQ and i/q are coordinates of the critical point. 
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a) 



b) 



c) 




FIG. 11: The global phase portraits for the phantom scalar field e = — 1 and the values of coupling constant: 

a) 1/6 < ^ < 3/16, b) 3/16 < ^ < 1/4, c) 1/4 < ^ < 3/10. In all cases there is present scenario of reaching 
the global attractor (a focus type critical point) from the unstable node. Note that in the case c) not all of 
the trajectories starting from an unstable node are reaching the de Sitter state, in contrast to cases a) and 

b) . 
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FIG. 12: The global phase portraits for the phantom scalar field e = — 1 and negative coupling constant 
f < 0: in figure b) f is greater, but still negative, than in figure a). It is easy to notice that in the limit 
^ 0^ we receive phase portrait from the Fig. [6^. for the phantom scalar field with minimal coupling. 



Using time transformation (jl3p we can calculate the scale factor growth along the trajectory 

(19) 



Alna= / V(o-)(l-e6^(l-60V'(o") )dCT, 
'o 



and the cosmological time growth 



f-af 1 />oo 1 

At= I — dlna = J — V'(ct) (1-66^(1 -60V'(f^)^)dCT. 



H 



(20) 



Linearized solutions are good approximations of the original system in the vicinity of the critical 
point. In what follows we assume that (V'(i) - ipo^ = [y{i) - Vof = - V'o) - Vo) = 0, 
a = 1 and yo > (see Fig. [T2]) . Then 



A In a = e(l - 6^) (V'(i) - V'o) V'o 



(21) 
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and 

I 

At = =el2^(l - 6^)(^/^(j) - V'o)V'o / V A exp {-yocr) + B exp {-2yQa) exp {-2yQ(j)da 

vrw^ Jo 

= ^.12^(1 - 60 (^(.) - ^o) { ( log ^ - 2 log ( + VaTB)) - 

- ^/B{A + B){-A + 2B){3A + 45)} (22) 

where ^ and S are positive constants 

A = -e{yo + 12^^-0) (2(1 - 30(V'(^) - V'o) + (y(^) - 2/o)) 

5 = £4(1 - GOiyo + 3CV'o)(V'{i) - V'o) 

For every case of existence of such a critical point (see Figs [2| [81191 1101 [TT] for unstable node and 
Fig. [12] for a stable node in the "physical region" ) growth of the scale factor and the cosmological 
time is finite. 

Now we calculate the first derivative of the Hubble function ([8]) with respect to the cosmological 
time at this point 

H = kn^y = ( r^-^ (23) 

2^ ^ Vl-e[^/;'2 + 12^V'V'' + 6^V^]/ ^ 
where a dot denotes differentiation with respect to the cosmological time and a prime with respect 
to the natural logarithm of the scale factor. Then after the elimination of second derivative of the 
field with respect to the natural logarithm of the scale factor we have 

1 - £6^(1 - eoV'^ (1 - £6^(1 - 60^2) (1 _ ^[^/2 + i2^^V' + eeV'^]) 

This expression at the critical point ipQ = ggg(i_gg-) ; V'o — ~^ i-6g ^ singular since the numerator 
is finite and the denominator is equal zero. 

The trajectories starting form this critical point corresponding to the singularities of finite scale 
factor seems to be most interesting. For such state appearing for both, canonical and phantom 
scalar fields we have curvature singularity because Hubble parameter is infinite but the scale factor 
assumes finite value. They are typical because the critical point is an unstable node (see Fig. [2] for 
an unstable node for the canonical scalar field and and Figs. [HI [U [10] and [11] for an unstable node 
and Fig. [12] for a stable node for the phantom scalar field) . 

Next we can proceed to the analysis of behaviour of system (I12p at the circle at infinity. In- 
troducing the polar variables in order to compactify the phase space by adjoining the circle at 
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infinity 



r r 

■0 = cos 9, y = sin ( 

1 — r 1 — r 



where r and 9 are the polar system coordinates, we receive the following dynamical system 

a— r(l — r)'^ = r(l — r) sin6'< cos6'(cos0 — sin 6*) ( —(1 — r)^ — e6^(l — 64)r^ cos^ 6* ) + 
d(T I Vk / 

+ e{l-Qi)r^ sin^ 9 cos 6* (sin 6* + 6^ cos 0)- 

'-(1 - rf - er"^ {{I- Qi)s\.n^ 9 + Qi{s\.n 9 + cos 9 f 

(e'^ ((1 - rf - e^r^ cos 6l(sin 9 + cos 9)) + 2r^ cos 6l(sin 6* + 6^ cos 61)) | , (25a) 

d^ I / 6 \ 

a— r(l — r)^ = cos 6'< — r^sin 0(sin 9 + cos9)(-{l- rf - e6^(l - 60^^ cos^ 9 ] + 
dcr I \K J 
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+ e(l — 6^)r^ sin^ 9 cos 0(sin + 6^ cos ( 
-( 

6 



[1-rf -er"^ ((1 - 6^) sin^ 9 + 6^ (sin 6* + cos 9f 
(^e^ ((1 - rf - EKir"^ cos 6l(sin 9 + cos 0)) + 2r2 cos 6'(sin 6* + 6^ cos 6*)) I . (25b) 



For identification of the critical points at infinity we can calculate the polar angle 9 from those 
equations simply by putting r = 1. Then we receive the following equation for the direction 9 
determining the localization of the critical points 

e3 cos^ 6* (2^ cos 6* + (1 - 4^ sin 9) ( sin^ + 6^ cos^ 9 + 12^ sin0cos0)=O. (26) 

We can notice that position of the critical points at infinity do not depend on the form of scalar field 
assumed, i.e. critical points are the same for both canonical and phantom scalar fields. In Table [III 
we have gathered positions and the character of these points. Note that there are specific values 
of the coupling constant ^ for which some critical points coincide. For example for minimal ^ = 
and conformal ^ = \ coupling critical points 2) and 3) in Table [TI] have the same location. This 
is the reason why in these cases the critical points at infinity are degenerated (i.e. the eigenvalues 
of linearization matrix calculated for the critical points located at infinity is identity zero). Such 
points are called multiple critical points. In our approach we treat the coupling constant ^ as a 
control parameter for which the bifurcation analysis is performed. The values of ^ different from 
minimal and conformal coupling split multiple critical points and remove the degeneration. 

Now we can simply calculate the value of the Hubble function ^ at critical points at infinity. 
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TABLE II: The critical points at infinity and their characters. 





Critical point 


existence 


eigenvalues 


1) 


cos = 


ve e R 


Ai = A2 = 


2) 


tanO — 


veeR\{^} 


Ai = £6^(1 - 60t^, A2 = -sUe^^ 


3) 


tan6' = -6^±v/6C(6^-l) 


^ < or ^ > i 


Ai = -el2e(l - 60(3(1 - 40 T 2^6^(6^ - 1)), 

A2 = -e6e(i - 60 (ee ± ^em - 1)) 



Using the polar coordinates with compactification with circle at infinity we receive 

-r^9i 9 cos^ 

H^\oo = -em^ . o„ , . , — — ^ . (27) 



Inserting the position angle of the critical points at infinity to the above expression we can conclude 
that at critical point 1) in Table [III = 0, for the second point = ~^'>tT''^ 2^(i-6i){3-i6^) ' 
finally for the point 3) the Hubble function is singular = 00. The most interesting critical point 
seems to be the second critical point because for the wide range of values of the parameter ^ the 
final state can be the de Sitter attractor with = const in spite of that the phase variables tp 
and ip' assume the infinite values. 



IV. CONCLUSIONS 



We study the dynamics of a scalar field with a simple quadratic potential function and non- 
minimally coupled to the gravity via ^Rip'^ term, where R is the Ricci scalar of the Robertson- 
Walker spacetime. We reduce dynamical problem to the autonomous dynamical system on the 
phase plane in the variables ^ and its derivative with respect to the natural logarithm of the scale 
factor. The constraint condition is solved in such a way that the final dynamical system is free from 
the constraint and is defined on the plane. We investigate the whole dynamics at the finite domain 
and at infinity. All the trajectories for all admissible initial conditions are classified, and critical 
points representing the asymptotic states (stationary solutions) are found. We explore generic 
evolutionary paths to find the stable de Sitter state as a global attractor and classify typical routes 
to this point. We study the effects of the canonical scalar field as well as the phantom scalar field. 
The following conclusions, as the results of our studies, can be drawn: 

1. The cosmological models with the quadratic potential function and non-minimal coupling 
term S^Rif)'^ are represented by a 2-dimensional autonomous dynamical system which is stud- 
ied in details on the phase plane (^/j, (see Table |T] for the critical points at the finite 
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domain). The shape of the physical region > does not depend on the form of the 
potential function, but only on the value of the coupling constant ^; 

2. We investigated the fixed points of the dynamical system and their stability to find the 
generic evolutional scenarios. We have shown the existence of a finite scale factor singular 
point (both future and past) where the Hubble function as well as its first cosmological time 
derivative diverge = oo and H = oo; 

3. For the phantom scalar field e = — 1 we found existence of a sink type critical point (i.e. 
a stable node or a focus, depending on the value of the parameter ^) which represents the 
de Sitter solution. Only for < 1/6 the evolutional paths avoid a past finite scale factor 
singularity (see for example Figs [6] and [8]) ; 

4. For the canonical scalar field e = 1 we found that for < ^ < 1/6 exists the critical point 
which corresponds to a past finite scale factor singularity for models with V{'ip) = i^m'^tlP' > 0. 
If < then this point corresponds to a future finite scale factor singularity (see Fig. [2]). 
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